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A fermion node is subset of fermionic configurations for which a real wave function vanishes 
due to the antisymmetry and the node divides the configurations space into compact nodal cells 
(domains). We analyze the properties of fermion nodes of fermionic ground state wave functions for a 
number of systems. For several models we demonstrate that noninteracting spin-polarized fermions 
in dimension two and higher have closed-shell ground state wave functions with the minimal two 
nodal cells for any system size and we formulate a theorem which sumarizes this result. The models 
include periodic fermion gas, fermions on the surface of a sphere, fermions in a box. We prove the 
same property for atomic states with up to 3d half-filled shells. Under rather general assumptions 
we then derive that the same is true for unpolarized systems with arbitrarily weak interactions using 
Bardeen-Cooper-Schrieffer (BCS) variational wave function. We further show that pair correlations 
included in the BCS wave function enable singlet pairs of particles to wind around the periodic box 
without crossing the node pointing towards the relationship of nodes to transport and many-body 
phases such as superconductivity. Finally, we point out that the arguments extend also to fermionic 
temperature dependent /imaginary-time density matrices. The results reveal fundamental properties 
of fermion nodal structures and provide new insights for accurate constructions of wave functions 
and density matrices in quantum and path integral Monte Carlo methods. 
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INTRODUCTION 

Let us consider a system of fermions described by a 
real wave function ^'(i?) where R denotes fermions co- 
ordinates R = (ri, r^r). Due to the antisymmetry, 
there exists a subset of fermion configurations for which 
the wave function vanishes and this subset is called the 
fermion node. The fermion node can be implicitly defined 
by ^'(i?) = 0, assuming that the nodal set does not in- 
clude configurations for which the wave function vanishes 
because of other reasons, eg, boundary conditions. In 
general, for TV spin-polarized fermions in a d-dimcnsional 
space, the fermion node is a {dN — 1)— dimensional mani- 
fold (hypersurface) . It is a well-known fact that for d > 1 
the antisymmetry alone does not specify fermion nodes 
completely. This is not difficult to understand since an- 
tisymmetry fixes only lower-dimensional coincidence hy- 
perplanes with dimensionalities {dN — d) . Therefore the 
fermion nodes and their properties are determined by in- 
teractions and many-body effects. 

The fermion nodes play an important role in quan- 
tum Monte Carlo (QMC) methods which belong to the 
most promising and productive approaches for studying 
quantum many-body systems. Let us briefly describe the 
basic idea of QMC and its relationship to fermion nodes. 
Consider a Hamiltonian H and a trial function \1/T(i?) 
which approximates the ground state ^'o(-R) of H within 
a given symmetry sector. It is straightforward to show 
that limT-^oo exp(— TiJ)\['T oc 4*0 where r is a real param- 
eter (imaginary time). This imaginary time projection 
can be conveniently carried out by simulating a stochastic 
process which maps onto the imaginary time Schrodinger 
equation. The wave function is represented by an en- 



semble of i?-space sampling points which are propagated 
according to G{R,R',t) = {R\exp{—TH)\R') and for 
large r the ensemble becomes distributed according to 
^I'o- Unfortunately, the straightforward application of 
this idea to fermionic systems encounters a fundamen- 
tal complication in the fermion sign problem (ll, |3i 3- 
The fermion sign problem makes QMC studies of large 
fermionic systems difficult since the statistical errors of 
fermionic expectation values grow exponentially in the 
projection time r and also in the number of particles N. 

One possibility for avoiding the fermion sign problem is 
to employ the fixed-node approximation which restricts 
the fermion node of the solution ^'o(i?) to be identi- 
cal to the fermion node of an appropriate trial func- 
tion ^fxiR)- The fixed- node approximation introduces 
an energy bias which scales as the square of the dif- 
ference between the exact and approximate nodes and 
therefore the accuracy of fermion nodes becomes very 
important. In particular, for the exact node one can 
obtain the exact energy in computational time propor- 
tional to a low-order polynomial in N. The fixed-node 
approximation has been very successful even for rather 
approximate nodes of commonly used trial wave func- 
tions which are based on Hartree-Fock (HF) determi- 
nants or a few-determinant post-HF expansions. The 
fixed-node QMC electronic structure calculations using 
HF or post-HF nodes usually recover between 90 and 95 
% of the correlation energy and energy differences agree 
with experiments typically within a few percent. Such 
encouraging results have been observed across many sys- 
tems such as atoms, molecules, clusters and soHds W\. 
It is remarkable that QMC methods have enabled us to 
" corner" the correlation energy problem into the last few 
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percent of the correlation energy using algorithms which 
polynomial scaling. 

Unfortunately, for many problems even a few percent 
of the correlation energy can be significant. Typical ex- 
amples are transition metal systems where the fixed-node 
error can be of the order of several eVs. It is therefore 
clear that better understanding of fermion nodes could 
be an important step forward for the QMC methodol- 
ogy and beyond. In addition, the fermion nodes are re- 
lated to other physical quantities and could shed light 
on other many-body phenomena which currently are not 
completely understood. 

In order to introduce the basic properties of the 
fermion nodes and to illustrate the problems involved we 
will first mention a few unique cases of interacting sys- 
tems for which the nodes are known exactly. These exam- 
ples include a few two- and three-electron atomic states, 
namely, triplets of He atom ^S{ls2s), ^P{2p'^) ^ and 
the exact node of a three-electron atomic state '^S{p'^) 
The wave functions of these high symmetry states 
can be parametrized by appropriate coordinate maps in 
which the node is described by a single variable. For ex- 
ample, in the case of He triplet ^S{ls2s) state the corre- 
sponding "nodal coordinate" is cos/3 — r+2 • I"i2/(?'i2''i2) 
where rfj = -I- rj, Vij = r.^ — rj. For the triplet ^P(2p^) 
the relevant variable is cosw' = Zq • (ri x assuming 
that the P state is oriented along the z-axis which is 
specified by the unit vector Zq 0, 0|. For the quar- 
tet state *<S'(p^) the node is captured by the variable 
cos a; = r 1 • (r2 X ra ) . The node in these systems is encoun- 
tered whenever the corresponding variable, cos /3, cos w' 
or cosw, vanishes. 

The nodal surface divides the space of fermion config- 
urations into nodal cells, sometimes also called domains. 
(The name nodal cell was introduced in the previous pa- 
per In algebraic geometry and topology nodal cells 
are usually called nodal domains, see, for example, Ref. 
01 . We will use both of these two expressions inter- 
changeably.) The few known examples of exact nodes 
mentioned in the preceding paragraph help to illustrate 
an important point, which has been conjectured for the 
fermionic ground states for some time pj, namely, that 
the ground state node divides the configuration space 
into the minimal number of two nodal cells: a "plus" cell 
with ^' > and a "minus" cell with \[' < 0. From the 
node equation cos/3 = of the ^5(1525) state we easily 
find that the "plus" and "minus" domains are given by 
the conditions ri > r2 and ri < r2, respectively. For 
the node given by cos lu' — the nodal domains are 
given by the orientation of zo,ri,r2: the three vectors 
are either left- or right-handed, corresponding to either 
"plus" or "minus" cell, respectively. Similarly, for the 
node cosw — the domains are given by the left- or 
right-handedness of the vectors ri,r2,r3. The minimal 
number of two nodal domains was found also for 2D and 
3D noninteracting spin-polarized homogeneous electron 



gas with periodic boundary conditions for up to 200 par- 
ticles using a numerical proof . 

Understanding the fermion nodes and their properties 
has become a challenge which might help to advance both 
practical calculations but also open a deeper insights into 
the properties of fermionic systems. One can envision two 
key problems: 

a) Topology of the fermion nodes, ie, how many nodal 
cells are actually present since this is of high importance 
for the fixed- node QMC approaches. Note that if the 
number of nodal cells is incorrect (typically higher than it 
should be since mean-fields such as HF have tendency to 
divide the configurations space into too many domains) 
then the QMC sampling around the artificial nodes will 
be very sparse. This could lead to large statistical fluctu- 
ations from poor sampling, and, possibly, to an effective 
non-ergodicity due to the finite-time projection time in 
practical calculations. 

b) Once the topology is correct, the accuracy of the 
manifold shape becomes important. This is an area where 
our insights are particularly limited since the exact nodes 
for large interacting systems are virtually unknown ex- 
cept for a few-particle special cases mentioned above. 

In our recent paper |l7l | , we have made some encourag- 
ing steps forward in trying to understand the topological 
issues and we have analytically derived a number of new 
results regarding the number of nodal domains. In par- 
ticular, we have shown that spin-polarized closed-shell 
ground states of noninteracting harmonic fermions in 
d — 2 and higher have the minimal number of two nodal 
domains. We have proved the same for spin-polarized 
atomic states with several electrons, both for noninteract- 
ing and HF wave functions. We have also shown that by 
imposing additional symmetries one can generate more 
than two nodal domains but that interactions can relax 
this "nodal degeneracy" to the minimal number of two 
domains such as in the case of the ^S'(ls2s3s) atomic 
state. 

For noninteracting spin-unpolarized systems, ie, with 
both spin channels occupied, the number of nodal cells 
is four since the wave function is a product of spin-up 
and spin-down Slater determinants (2x2 = 4). (Here 
and later on we assume that the Hamiltonian does not 
include spin-dependent terms so that the particle spins 
are conserved. The wave function is then a product of 
two determinants which depend only on the spatial de- 
grees of freedom |2j.) In the last few years, studies of a 
few-particle systems have revealed that interactions and 
many-body effects do affect the nodal topologies and can 
change the number of nodal cells 0, 0. In particular, 
for the case of Be atom it has been found that the nonin- 
teracting/HF four nodal cells of the singlet ground state 
change to the minimal number of two due to the elec- 
tron correlation |^ and qualitatively the same has been 
observed in other systems ou'' recent paper 0| 

we have found that this is a rather generic property of 
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ground states in interacting systems. With some condi- 
tions, we have explicitly demonstrated that interactions 
lift the "nodal cell degeneracy" in spin-unpolarized sys- 
tems and smooth out the four nonintcracting cells into 
the minimal two. We have shown this for 2D harmonic 
fermions in closed-shell states of arbitrary size using a 
variational Bardeen-Cooper-Schrieffer wave function. 

In this work we further advance these ideas. We an- 
alyze the fermion nodes in a number of other paradig- 
matic fermionic models: homogeneous electron gas with 
periodic boundary conditions, particles in an infinite 
well and on the surface of a sphere. For all these sys- 
tems we prove that in the spin-polarized nonintcract- 
ing closed-shell ground state the fermion nodes divide 
the configuration space into minimal two cells for ar- 
bitrary number of particles. We also extend our pre- 
vious proof for spin-unpolarized systems demonstrat- 
ing that interactions smooth out multiple nodal cells of 
noninteracting/mean-field wave functions into the mini- 
mal two for more systems such as homogeneous electron 
gas and 3D harmonic oscillator. These results contribute 
to our understanding of the fermion nodes and their im- 
pacts on the accuracy of wave functions with direct im- 
plications for the QMC methods. 

In the last sections we show how in periodic spin- 
unpolarized interacting fermion gas the pairs of particles 
can wind around the box without crossing the node what 
points towards the connection of nodes to the transport 
and to the existence of other quantum phases such as su- 
perconductivity. Finally, we then generalize the results to 
the temperature density matrices with the implications 
for path integral Monte Carlo methods . 

II. GENERAL PROPERTIES OF FERMION 
NODES. 

Let us introduce the basic properties of fermion nodes 
as they were studied by Ceperley some time ago. 

a) Nondegenerate ground states wave functions fulfill 
the so-called tiling property. Let us define a nodal cell 
r2(i?o) as a subset of configurations which can be reached 
from the point Rq by a continuous path without crossing 
the node. The tiling property says that by applying all 
possible particle permutations to an arbitrary nodal cell 
one covers the complete configuration space. Note that 
this does not specify how many nodal cells are there. 

b) If m nodal surfaces cross then the angle of crossing 
is tt/to. Furthermore, the symmetry of the node is the 
same as the symmetry of the state. 

c) It is possible to show that there are only two nodal 
cells using the following argument based on triple ex- 
changes. Let us first introduce the notion of particles 
connected by triple exchanges. We will call the three 
particles i, j, k connected if there exists a triple exchange 
path i j, j ^ k,k ^ i which does not cross the node. 



If more than three particles are connected then they can 
form a connected cluster. An example of six particles 
connected into a single cluster is sketched as follows:«$«{. 
If there exists a point i?j such that all the particles are 
connected into a single cluster then ^'(i?) has only two 
nodal cells. This can be better understood once we real- 
ize the following two facts. First, any triple permutation 
can be written as two pair permutations. Therefore the 
connected cluster of triple permutations enables to realize 
any even parity permutation without crossing the node. 
That exhausts all permutations which are available for 
cell of one sign since the wave function is invariant to even 
parity permutations. Second, the tiling property implies 
that once the particles are connected for Rt the same is 
true for the entire cell. By symmetry, the same argu- 
ments apply to the complementary "minus" cell which 
correspond to the odd permutations. More details on 
this property can be found in the original Ref. [^. 

III. NONINTERACTING SPIN-POLARIZED 
FERMIONS. 

III. a. Homogeneous electron gas. 

We consider a system of spin-polarized noninteracting 
fermions in a periodic box in d dimensions. The spatial 
coordinates are rescaled by the box size so that we can use 
dimensionless variables and the box becomes (— 7r,7r)''. 

We first analyze the fermion nodes for d = 1 since 
the result will be useful in subsequent derivations. We 
consider a system with N = (2kp + 1) particles. In our 
ID dimensionless units the Fermi momentum becomes 
an integer, kp = 1,2.... The one-particle occupied states 
are written as (f>n{x) = e™^, n = 0, ±1, ...,zLkp and the 
spin-polarizcd ground state is given by 

^lDil,-,N)^det[(j)n{Xj)] (1) 

where xj is the j-th particle coordinate and j — 1, TV. 
We factorize the term exp(— j/c^ Xj) out of the deter- 
minant so that the Slater matrix elements become pow- 
ers of Zj = e*^J . The resulting Vandermonde determinant 
can be written in a closed form and after some rearrange- 
ments we find 

iV) = e-'"- Y[iz, ~ Zk) = 

j>k 

= A^o n sin(2;jfe/2) (2) 
j>k 

where Xjk — xj — Xk and /io is a constant prefactor which 
is unimportant for our purposes. (In the derivations be- 
low we will be using the letter jj, for denoting prefac- 
tors which are either constants or nonnegative functions 
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in particle coordinates and therefore do not affect the 
nodes) . The derived wave function has the following im- 
portant properties. 

a) The fermion nodes appear at the particle coinci- 
dence points. This implies a well-known result that the 
ground state wave function in ID have A^! nodal cells 
since any particle permutation requires crossing the node 
at least once. In addition, this also means that any 
fermion configuration which preserves the particle order 
is contained within the same nodal domain/cell. 

b) The wave function is invariant to translations and 
cyclic exchanges of particles. Due to the periodic bound- 
ary conditions this includes also winding the system 
around the periodic box, ie, the translation of all par- 
ticles by 2tt. We can formalize this by introducing an 
operator which translates all the particles along the 
cc-axis as Xj Xj + a so that the translation invariance 
can be written as 



where kp in 2D is not necessarily an integer. Similarly to 
our previous paper ^3 1 '^^ show that the spin-polarized 
electron gas for closed-shell ground states have only two 
nodal cells. 




(a) (b) (c) 

FIG. 1: Positions of five fermions in the 2D periodic box 
(— 7r,7r)'^. (a) Particles aligned horizontally, along the lines 
y = Co,CiiC2- (b) Particles aligned vertically, along the lines 
X — ^o,Ci)C2. (c) Particles aligned in both directions, posi- 
tioned on a lattice. 



t:^M^,...,n)^^M^,...,n) (3) 

Similarly, the wave function is invariant to cyclic ex- 
change of all the particles given hyj^j + l,j — 1, ...,N 
and iV + 1 1. This action is carried out by C^i op- 
erator where the notation means the exchange by one 
site is in the -|-a;-direction. Clearly, the inverse operator 
[C*!.]^]"^ = C^i- Note that here we have assumed that 
N is odd, in agreement with our definition. The invari- 
ance holds only for TV odd since then the cyclic exchange 
is equivalent to an even number of pair exchanges. For 
the sake of completeness, we consider also TV even, when 
the cyclic exchange can be replaced by an odd number of 
pair exchanges, resulting in the wave function sign flip. 
In general, for the cyclic exchanges we can therefore write 

C^.^Ml, TV) = i^lf+'^Mh N) (4) 

c) Assuming the particle positions are all distinct, the 
cyclic exchange path can be chosen in such a way that it 
does not cross the node. This is easy to accomplish by 
maintaining finite distances between the particles along 
the path trajectory. Let us parametrize the exchange 
path by a parameter t so that the path starts at i = 0, 
the path is completed at t = 1 and the exchange path 
operator is then denoted as C':^i{t) with < t < 1. The 
fact that the path does not cross the node (ie, the path is 
contained within the same nodal cell) can be then written 
as 

|C^i(i)*iD(l,...,TV)| >0, 0<t< 1 (5) 

where, of course, TV is assumed to be odd. 

Let us now derive the wave functions and generalize the 
results for fermions in 2D. The one-particle states in 2D 
are (j)nm{x, y) = e*("^+™i'). The states are occupied up to 
the Fermi momentum kp so that we have + rn? < kp, 



The proof is by induction, therefore let us first consider 
kp ~ 1 with five particles occupying {1, 6"*=", e^'^} one- 
particle states. We place the particles as in FigQJi so that 
the coordinates are given by ri = (a:i,Ci), r2 = (a:2,Ci)j 
i"3 = (a;3,Ci), i"4 = (0:3,(2), rs = (0:5,(0 )■ By eliminating 
the terms with (1 the wave function can be factorized as 
follows 
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= ^i'o^lD{l, 2, 3) sin(Cio/2) sin(C2o/2) sin(C2i/2) (6) 

where /xo,/Lto are irrelevant constant prefactors. We have 
obtained a product of ID wave function and terms with 
distances between the lines y = Co, Ci, C2- Note that anal- 
ogous result can be found for the configuration sketched 
in Fig. lb so the particles are aligned in parallel to the 
y-axis with the wave function given by 

*2d(1,--,5) = /^o*id(5,2,4)x 



X sin(fio/2) sin(6o/2) sin(f2i/2) (7) 

If the particles are aligned in both directions as in Fig. Ic 
both expressions apply. From the derived wave function 
it is clear that the node is encountered when particles 
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lying on the same line reorder or when the lines of parti- 
cles cross each other, eg, = Co- The wave functions for 
the configurations outhned in Fig^possess an important 
property. Consider in Eqs. IHlElwith groups of 

particles positioned on the corresponding lines. Assum- 
ing the group of particles is allowed to move only along 
the given line, one can consider this to be an effective 
ID subspace. The ID wave function for such a subspace 
obeys all three conditions derived for ID case, ie, Eqs. 13 
0] and El For example, for ^'id(1,2, 3) which appears in 
Eq. El we have |C^i(t)4'ir,(l, 2, 3)| > 0. This is easy to 
check also analytically. We fix the coordinates in Fig. 1(a) 
as follows xi = — 27r/3,a;2 = 0,2:3 = 27r/3 and then we 
can carry out a " synchronized" cyclic exchange using the 
translation by 27r/3 in x-direction. The wave function 
is constant along the translation/cyclic exchange path so 
that we have 

C^i(0«'id(1,2,3) = *ii5(1,2,3), 0<t<l (8) 

because of the translational invariance. 

Finally, we are ready to show that there are only two 
nodal cells in this five-particle ground state. Assume that 
the particles are positioned as in Fig.lc. Consider the 
following sequence of four exchanges, C^iC^iC^iC^i, 
where the operators act from the right, ie, C^^, acts first. 
We denote this exchange symbolically as 



/~ix /~iy /~ix /~iy 
'-'+1^+1^-1^-1 



(9) 



The exchanges are performed only for the particles lying 
on the lines Ci;Ci in corresponding directions, eg, C^j^ 
acts on particles along y = Ci while cycles particles 
along X = £^1 . It is easy to check that this results in a 
triple exchange 423 342 while the particles 1 and 5 end 
up in their original positions. It is therefore clear that 
particles 2,3,4 are connected by this triple exchange. It is 
straightforward to show that C''tiC^_iC^iC^i exchanges 
particles 1,2,5 so that we can conclude that all five parti- 
cles are connected by triple exchanges into a single cluster 
and there are only two nodal domains. 

Now we need to extend the arguments and find the 
wave functions for general cases with occupied states 
within the Fermi disk with an arbitrary kp. It is in- 
structive to first derive the wave functions for kp = \/2 
with 9 particles occupying states {1, e^", e^'^, e^'^^^^^^} 
and for kp — 2 with 13 particles with additional states 
|g±i2i:^ g±42j/|_ fpj^jg enable us to understand how to 
write down the wave function for a general case while 
avoiding rather tedious notations which would appear in 
the direct derivations. 

Let us position the particles into a pattern which mim- 
ics the lattice of the occupied wave vectors in the recipro- 
cal space, see Figj^l Using appropriate algebraic arrange- 
ments we subsequently factorize the lines of particles as 
given by 
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FIG. 2: Positions of fermions in the 2D periodic box in real 
space, (a) Positions of the nine-fermion system, (b) Positions 
of the 13-fermion system. See the text for details. 



= Mo*id(1, 2, 3)*2d(4, 9) sin3(C2i/2) sin3(Cio/2) 



/i[,*iD(l, 2, 3)*id(4, 5, 6)*iz)(7, 8, 9) x 



X sin3(C2i/2) sin3(Cio/2) sin3(C2o/2) (10) 

where /iOiMo constants. Similarly, for kp — 2 we can 
first factorize the line with the largest number of particles 
on y = (^2 (see FigjSJ 

^2d{1, ■-. 13) = ^^0^-1^(1, 5) sin(Co2/2) sin3(Ci2/2)x 



X sin3(C32/2) sin(C42/2)*2i5(6, 13) 

and then factorize out the particles lying on another line 
2/ =: Cs to obtain 



^'2Z3(1,-,13) =^^0 



nsin"^(C,2/2) 

j¥2 



n sin"^(C,i/2) 



*ii5(l, 5)«'iz3(6, 7,8)^-2^(9, 13) (11) 

where rij is the number of particles lying on the line 
y = Q. We have obtained a product of ID wave func- 
tions, terms with distances between the lines and the 2D 
wave function with lower number of particles positioned 
in the same type of pattern is in Fig^ Obviously, we 
can further factorize \E'2_d(9, 13) using EqjSl until we 
end up with ID factors only. Using these insights into 
the recursive form of the wave function, for a general case 
with AI + 1 lines we can write 



*2d(1,...,9) = 
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where Ih 



. , in'u denotes the labels of particles ly- 
ing on the line y = Ck- In addition, analogous expression 
can be found if we factorize along the x ~ lines, the 
only difference being replacement by and corre- 
sponding replacement of particle sets in '^m wave func- 
tions. 
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FIG. 3: Illustrations of the position patterns in 2D periodic 
box for the size increases from to fc^ -I- Akp- The parti- 
cle layouts mimic the occupied states in the reciprocal space, 
however, the particles are positioned in the real space. The 
additional particles are in grey/yellow, (a) The additional 
particles positioned along x— and j/— axis directions, (b) The 
additional particles positioned along diagonals. Lines x = 
and y = C,i are used for the proof that the particles a, b, c are 
connected by a triple exchange. See the text for details. 

We are now ready for the induction step. Consider a 
spin-polarized closed-shell ground state with a given kp. 
For this wave function we assume that all the particles 
are connected by the triple exchanges, ie, there are only 
two nodal domains. Let us increase kp kp + Akp until 
the Fermi disk includes the next unoccupied star of states 
4>nm{x,y) for which /ci? < {ji^ +m?Y/'^ < kp + Akp. The 
system size increases by the corresponding number of ad- 
ditional particles. Assuming the particles are positioned 
in the real space in the same pattern as the occupied k- 
points in reciprocal space (Fig|31) , the additional particles 
will appear at the borderline of the disk in the real space. 

The two basic possibilities how the additional particles 
are positioned are given in the Fig|31 We need to show 
that these additional particles are connected to the origi- 
nal particles by the triple exchanges. This can be demon- 
strated by the sequence of the four cyclic exchanges which 
we have used above for the five particle case. It involves 
particles on the lines x — S,k and y — Q schemati- 
cally drawn in Fig|2| Consider the exchange C^|«_j^ = 
C^iC^iC^iC^i where the cyclic exchanges in x direc- 
tion are applied only to particles on the line x = and, 
similarly, the cyclic exchanges in y direction are applied 
only to particles along y = Q. Note the wave function is 



invariant to cyclic exchanges since the number of particles 
along each line is odd for any closed shell state. It is then 
straightforward to find out that C^^^i exchanges parti- 
cles a, 6, c while the rest of the particles remains intact. 
Similar exchanges can be carried out for all additional 
particles. Finally, this shows that the additional parti- 
cles are connected to the particles of the wave function 
with size kp and finalizes the proof. 

The proof can be extended into 3-D and higher dimen- 
sion by positioning the particles onto an appropriate 3-D 
pattern which in real space mimics the occupied states in 
the Fermi sphere in the reciprocal space. This is possible 
due to the fact that with proper positioning of particles 
one can subsequently factorize the Slater determinant 
along hyperplanes, planes and lines. Using arguments 
similar to the 2D case we can perfrom cyclic exchanges 
without crossing the node in 3-D and higher dimensions. 
Therefore the proof for higher dimensions is essentially 
the same. 

In many cases, the proof can be extended to open 
shells. If an open-shell state is degenerate it is neces- 
sary to fix the ambiguity in the nodes, for example, by 
considering wave functions which transform according to 
an appropriate symmetry subgroup 0, In general, 
however, depending on symmetries, the number of de- 
generate states, etc, one cannot rule out possibilities of 
states with the number of cells beyond the minimal two 
(this will be investigated in the next paper). 

This concludes the arguments and the proof that for 
d > 1 the noninteracting spin-polarized fermion gas 
closed-shell ground states in periodic boundary condi- 
tions have only two nodal cells. 



III. a. Fermions on the surface of a sphere. 

Spin-polarized free fermions on the surface of a sphere 
are described by a Slater determinant with one-particle 
states being the spherical harmonics Yim.{'d,(p). The 
spherical harmonics are polynomials in variables cos?? 
and sini^e**'''. The Slater matrix elements can be re- 
arranged to monomials in these variables and it is then 
straightforward to factorize the determinant into similar 
form as obtained for the homogeneous fermion gas or for 
the harmonic oscillator • Let us assume that the par- 
ticles are positioned as sketched in Fig.^Ji. Using familiar 
expressions for the first few spherical harmonics we get 
for the ^S{sp^) state 

^!,p3{l,...,A) ^ fioiui-uo)vl Yl sin(v5y/2) (13) 

2<i<j<4 

where we have denoted u — cos 1} and v — sin -d and fxo is 
a constant. Similarly, for the S{sp^d^) state with nine 
fermions we obtain 



(1,...,9) = Mo*sp3(l,...,4)x 
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X (m2 - uo){u2 - ui)''f;^*iD(5, 9) 



(14) 



where we have denoted N) — Yij^k siii((pjfe/2) 

in agreement with Eq. [5] 



*3 



♦7 



♦7 



*6 



"0 

(a) 



(b) 



FIG. 4: Positions of the particles on the sphere surface for 
the state S{sp'^cf') using coordinate system u = cos(i?) and 
< f < 2n. The particles are positioned along the lines with 
constant u = Uj. For larger systems the patterns are the same 
with particles corresponding to increasing angular momentum 
/ lying on the lines with increasing ui. (a) Particle positions 
used for derivation of the wave function factorization, (b) 
Particles after appropriate shifts along ip and a subsequent 
rotation. 



For an arbitrary size closed-shell 5 symmetry state the 
one-particle states are occupied up to the angular mo- 
mentum L with the corresponding number of particles 
N = I]f=o(2^ + 1) ^ {L + If. Using the previous two 
examples we can directly write the factorized wave func- 
tion assuming the particle positions follow the pattern 
in Fig. We first factorize the line u — ul and then 
recursively the rest of the lines so that we can write 



4'(l,...,iV) = //0*2D(l,...,Af//L)x 



n 

0<i<L 



L-1 



fj.{vi, ...,vl) Yi 



k=0 



j<L-k 



(15) 

where Ik is the subset of particles lying on the line u — u^. 
In the final expression the powers of Vj , such as w^^^^^'' in 
the preceding line, are absorbed into zjl). This 

prefactor is nonnegative and does not affect the nodes or 
the wave function rotational invariance. 

We use the properties of ID wave functions and the 
rotational invariance to build the proof of the two nodal 
domains. First, note that for four particles in ^S{sp^) 
state it is easy to show [ivj that there are only two do- 
mains. A little bit of algebra shows that the node is 



encountered when all four particles lie on a circle result- 
ing from a plane cutting the sphere. Clearly, it is easy 
to position the four particles on the sphere so that the 
triple exchanges do not cross this node. For the induc- 
tion step assume that for the size L with N = {L + i f 
the particles are connected. We need to show that for 
the size with L ^ L + 1 with additional 2{L + 1) + 1 
particles, the additional particles are connected as well. 
Using the fact that particles can be shifted along the ip 
coordinate and rotated, one can reposition the particles 
as illustrated on the ^°S{sp^d^) state, Fig. By ap- 
plying the factorization to this particle arrangement we 
see that the additional particles are connected and the 
argument applies to an arbitrary size closed-shell state. 



III.c. Fermions in a box. 

Let us assume a system of fermions in a box (0, ir f 
with the condition that the wave function vanishes at 
the boundaries. For d — 1 the one-particle states are 
4'n{x) = sin(na;), n = 1, 2, .... Note that this can be writ- 
ten as 4>n{x) = sin(a::)J7„_i(cosx) where Un is the n-th 
degree Chebyshev polynomial of the second kind. We 
map the variables Xi to pi = cos(xi) so that pi G (—1, 1) 
for I = 1, iV. Note that the map x — > p is a homeomor- 
phism (ie, it is bijective and continuous with its inverse). 
Homcomorphisms preserve topologies, eg, the ordering of 
points, so that Xa < Xb <^ pa < Pb- Using this map we 
can write the wave function for N fermions in the \D 
box directly as 



'J'iD(l,...,Af) =det[0„(a;O] 



^0 J|sin(a;fe) \\{pi - Pj) 

k i>j 

(16) 

where fiQ is a constant. 

In the 2D box the one-particle states are given 
by (l)nm{x,y) = sin{x)sm{y)U„-i{p)Um-i{q), n,m = 
1,2,... where we denoted q = cos{y). If we absorb 
sin(x) sin(i;) into the common prefactor /i, the Slater ma- 
trix elements become monomials of the type p"'~^q™~^. 
The states which are occupied lie within the quarter of 
the Fermi disk -|- < kp and n, m > 0. Assume that 
the particles lie on a lattice in the space of variables p, q 
so that they are positioned on M horizontal and vertical 
lines (see, eg, the upper right quadrant in Fig. 3c). Using 
the techniques presented above and in the previous paper 
|l7| we can write down the wave function as 

^-2^(1,...,^^) = 



A/-1 



= Kxi,yi, ■■■,XN,yN) Yl 



k=l 



^lD{Ik)\{{pk-PiT 

j>k 



(17) 
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where the prefactor is a nonnegative function. Similar 
expression can be written down by factorizing along the 
lines q ^ Qi- The wave function therefore factorizes in a 
manner very similar to the harmonic oscillator and peri- 
odic fermion gas. Therefore the proof of the two nodal 
domains above can be constructed in a similar manner. 
First, it is not too difficult to show that the there are only 
two domains in the system with three particles. Next, 
one assumes that this is the case for a system with M 
lines and then it is straightforward to show that the same 
applies to the system with M + 1 lines. Note that EalTTI 
suggest that apparently the wave function has both ro- 
tation and translation invariance. This is strictly not 
true, due to the nonlinearity of the map {p,q) {x,y). 
What is true, however, is that node is not crossed during 
rotations and translations since (x, y) {p, q) is home- 
omorphic. Therefore rotations and translations change 
the wave function values but not the sign since there is 
no reordering of either the lines or the particles along the 
lines. That enables us to use rotations and translations 
to prove the connectedness assuming that we take addi- 
tional care to keep the paths contained within the box. 
The proof then follows similar line of arguments as pre- 
sented before. It is easy to show that the three-particle 
state has only two nodal cells. Place the particles on 
a lattice and position the system into the center of the 
box. We assume that the lattice constant is sufficiently 
small so that translations by one lattice constant would 
not push the particles out of the box. It is then straight- 
forward to show that the four translations T^^T^^T^Ty 
exchange three particles at the edge of the system and 
by similar exchanges one can show connectedness of all 
particles for arbitrary size. Therefore the ground state 
closed-shell wave functions for particles in a box have the 
same generic properties as the fermionic models studied 
in preceding sections. 



IV. MINIMAL NUMBER OF NODAL CELLS 
THEOREM FOR SPIN-POLARIZED 
NONINTERACTING SYSTEMS. 

Using the results and proofs derived in this section and 
also in the previous paper [l^ we can write the following 
theorem. 

Theorem. Consider noninteracting or mean-field spin- 
polarized fermions in d > 1 with an exact wave function 
given by a Slater determinant of one-particle states times 
a prefactor which does affect the fermion nodes. Let the 
one-particle states are such that the Slater matrix ele- 
ments can be rearranged to monomials in particle coor- 
dinates or in coordinates transformed by a homeomor- 
phic map. Then for an arbitrary size closed-shell ground 
state the corresponding wave function has the minimal 
number of two nodal domains. 

The theorem covers a number of paradigmatic models 



and it is quite suggestive to think about this as being 
related to general mathematical properties of zeros of 
functions defined through determinants. In fact, the fac- 
torizations which enabled us to prove the two nodal cells 
property is directly related to the properties of multiple 
hyperplane configurations and to the multi-variate Van- 
dermonde determinant theorem which can be found in 
mathematical literature in various contexts 0, 0| . 

Considering that we restricted the proofs to noninter- 
acting systems and we relied on the fact that the ma- 
trix elements are monomials, it is useful to consider cases 
which go beyond such a framework. This line of thought 
leads us to the following interesting questions: 

i) Is the two nodal cell property valid for noninteracting 
cases with Slater matrix elements not reducible to mono- 
mials ? 

ii) What is the impact of interactions ? 

A tentative answer to the first question is given in 
the next section where we show that atomic states in 
Coulomb potential, which cannot be reduced to mono- 
mials due to the shell structure, exhibit the same prop- 
erty. We will not investigate the impact of interactions 
for spin-polarized systems in this paper. However, we 
will study and prove the two nodal cells for perhaps even 
more important cases of interacting spin-unpolarized sys- 
tems in the section VI. 

V. SPIN-POLARIZED ATOMIC STATES. 

In the previous paper we have proved that the 
atomic spin-polarized state Is2s2p^ has two nodal cells 
for both noninteracting and HF wave functions. The 
proofs for atomic states are more involved since for 
Coulomb potential it is more difhcult to find appropriate 
factorizations. The main complication is that orbitals 
in different subshells and angular momentum channels 
have, in general, different radial dependences which can- 
not be all factorized out of the determinant into a com- 
mon prefactor. Nevertheless, it is possible to demon- 
strate the two nodal cells for atomic states for several 
spin-polarized (half-filled) main subshells (and possibly, 
for all the states relevant for the periodic table of ele- 
ments). We will illustrate the idea of the proof on the 
spin-polarized ^'^S'(ls2s2p^3s3p^3d^) state with 14 elec- 
trons and then point out how the proof can be extended 
to larger systems. The onc-particle orbitals are pis(r), 
P28{r), P2p{r)x, P2p{r)y, P2p{r)z,...,p3d{r){2z'^ - x'^ - y'^), 
etc, and we use dimensionless coordinates which are 
rescaled as r ^ Zr/ao with Z being the nuclear charge 
and flo the Bohr radius. The wave function is given by 

*at(l,...,14) = 



det{pj,, y, z, pI^, (j)^^^, ^^^^s , ...} (18) 
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where pl^{r) = pis{r)/ p2p{r), p^^(r) = P2s{r) / P2p{r) and 
(l^lpxi^) = xp3p{r)/p2p{r), (f>laz2{r) = p3dir){2z^ - 

— y'^)/ P2p[r), etc, where we factorized the nonneg- 
ative radial function P2p{r) out of the determinant. 

We wiU show the connectedness of all the particles in 
two steps. We first demonstrate that the Slater determi- 
nant can be factorized into subdeterminants correspond- 
ing to subshells Is, 2s2p^ and 3s3p^3d^. That will enable 
us to show that the particles within each subshell are con- 
nected. In the second step we show, that the particles 
can exchange between the subshells. For this purpose we 
will specify the explicit paths and carry out a straight- 
forward numerical check that there is no node crossing 
by tracing the wave function along the paths. 

In order to show the factorization into shells we posi- 
tion the particles as follows: particle 1 is in the origin, 
particles 2 to 5 are on the surface of a sphere with the 
radius rja and particles 6 to 14 are are on the surface of a 
sphere with the radius rji,. The radii rja and ij}, are given 
by the radial nodes of the orbitals P2s(i') and P3p(r), ie, 
PisiVa) = and pzp{rih) = 0. The Slater determinant can 
be written in the form 



^'at(l,...,14) =det 



A D G 

B E' E 
C F H 



(19) 



where the block matrices A to H are given as follows 



A = 



pL(0) 
pL(0) 











ai 




ai 








X2 2^3 2^4 



2/2 

Z2 



2/3 



2/4 



2/5 



(20) 



B 



P3s(0) 







D 



03 03 


03 03 




'/'3p.(2) - 


-^3^.(5) 


03p,(2) - 


<Plpy{^) 


03p.(2) 


r^pzi^) 


L.^(2) - 






L.42) ... 






3(i£cy(2) ... 


4>lcixyi^) 




Myzi'^) ... 


r^yzi^) 






^3dxzi5) _ 





(21) 



(22) 



bi 


bi 


bi 


bi 


62 


62 


b2 


b2 


xe 






Xg 


2/6 






2/9 


Z6 






Z9 



(23) 



E : 



bs bs 63 &3 63 









(24) 



bldxyi^) 
Kdyzi^) 
^Idxzi^) 



r^dx^i^) 

<t>ldxy{^) 
r^dyzi^) 
r,dxzi^) 



(25) 



G = 



bi 


bi 


61 


bi 


bi 


b2 


62 


&2 


b2 


b2 


xio 








Xi4 


2/10 








2/14 


2:10 








Zi4 



(26) 



H 



Pidz^-m 
^sdx^m 

^idxyim 
r^dyzm 

Kdxzm 



^3..^ (14) 
^3d.^(14) 
*3d.,(14) 
^3d,.(14) 

^Idxzi^^) 



(27) 



and ai = plsiVa), 03 = P^iVa), &i = PuiVb), &2 = 
/92s (?7f)), ^3 — PssiVb)- The block matrix E' is the same 
as E except that it has only four columns instead of five. 

The following three row additions in the Slater matrix 
given by Eg 1191 allow for factorization into subdetermi- 
nants. We first eliminate elements ai in the first row by 
adding an appropriate multiple of the fifth row. Simi- 
larly, we eliminate the terms with bi in the first row by 
adding a multiple of the second row. This leads to the 
first row having only one nonzero element A'n = plg{0) 
— aiP3s(0)/a3 —(61 — ai53/a3)p2s(0)/^2 and reduces the 
Slater matrix by the first row and the first column. Fi- 
nally, by adding a multiple of the second row to the fifth 
row, both E' and E become zero matrices and we can 
write 



«'at(l,...,14) = A'lidetB'det 



D' G' 
F H 



(28) 



where B' is the matrix B without the first column, while 
D', G' are the matrices D, G without the first row. The 
three factors in the expression above correspond to the 
three main subshells namely, Is, 2s2p'^ and 3s3p^3c?^ with 
the dependence on positions of particles 1, 2-5 and 6-14, 
respectively. The key point is that each of the subshells 
represents a system of particles on a sphere in the S*- 
symmetry state and for such cases we proved that the 
particles are connected. That concludes the argument 
that the particles within each subshell are connected. 

More difficult part of the proof is to show that one can 
exchange the particles between the subshells. We were 
able to prove analytically the exchange between Is and 
2s2p'^ in our previous paper. For larger cases the ana- 
lytic proof becomes very tedious and it is much more ef- 
ficient to evaluate numerically the determinant along the 
following exchange paths. The particles are positioned 
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FIG. 5: Positions of 14 electrons for the proof that triple 
exchanges connect the particles in the S{ls2s2p^3s3p'^3(f') 
state. Particles 2 to 5 lie on the spherical surface with the 
radius equal to the radial node of p2s{r) orbital. Particles 6 
to 14 lie on the spherical surface with the radius equal to the 
radial node of p2p{r). 



as sketched in Fig|Slwith the coordinates given by ri — 
(0,0,0), V2 = (77a, 0,0), r2 = (0,77a,0), r4,5 = (0,0, ±77,), 
rg - (776/\/2,77fc/\/2,0), rr.s = (±??6/\/2, T776/\/2, 0), 
rg.io = (0,0, ±?7f,), = (77b/\/2, 0, ±?7b/\/2) and 

ri3,i4 — (0, 77b/\/2, ±?7b/\/2). There exists enough of 
triple exchanges of neighbouring particles, along the sides 
of corresponding triangles, to connect all the particles 
into a single cluster. For example, the wave function 
values for exchanges 123 231, 236 362 (and other 
exchanges for illustration) are given in Fig. El 




FIG. 6; Wave function of the i'^S'(ls2s2p^3s3p^3d'^) state 
(arb. u.) along triple exchange paths parametrized by t. Par- 
ticles are positioned as given in Fig.|3 The plotted exchanges: 
1,2,3 (fuU line); 2,3,4 (dashed line); 2,3,6 (dotted line); 6,7,9 
(dashed-dotted line); 9,11,13 (double dashed-dotted line). By 
symmetry, these non-crossing exchanges connect all the par- 
ticles of this state. 

This is enough to show that the three subshells are con- 
nected since many other exchanges are identical due to 
the symmetries in particle positions and the wave func- 
tion S symmetry. For the illustration in Fig. we used 
the noninteracting radial orbitals. For HF orbitals one 
gets the same qualitative picture since the basic spatial 



properties of the noninteracting and HF orbitals are qual- 
itatively the same (ie, ordering of radial nodes of pni (r) 
orbitals, behaviour at nucleus, tails, etc). The proof uses 
only the fact that some of the radial nodes of the one- 
particle orbitals are ordered as in the non-interacting case 
so it can be extended to HF wave functions as well. 

One can expand the proof to larger systems, such as for 
the state with occupied fourth main subshell 4s4p'^4ci^4/'' 
and beyond. The factorization is similar: the particles in 
the fourth subshell are positioned on the spherical sur- 
face with the radius equal to the radial node of P4d(r) 
orbital. Somewhat long but straightforward algebraic re- 
arrangements show that the Slater determinant of the 
30 X 30 matrix can be reduced to product of subdetermi- 
nants corresponding to the subshells. For the purpose of 
this paper we do not deem necessary to go through the 
explicit proof. 

VI. INTERACTING SPIN-UNPOLARIZED 
FERMIONS. 

It is straightforward to understand the fermion nodes 
of noninteracting spin-unpolarized or partially polar- 
ized systems with more than one-electron in each spin- 
channel. The wave function is a product of spin-up and 
-down Slater determinants and the number of nodal cells 
is the product of the number of cells in each subspace. 
(As mentioned before, we assume that the Hamiltonian, 
even with interactions considered below, does not include 
terms with spin so that particles can be assigned to the 
spin subspaces.) For the ground states with two nodal 
cells in each subspace we get 2x2 = 4 nodal cells. This 
is the nodal structure of a single configuration Hartree- 
Fock wave function and as such has been used in many 
fixed-node QMC calculations. However, analyses of small 
interacting systems revealed that this is not correct and 
interactions do change the nodal topologies and the num- 
ber of nodal cells. For the first time this has been demon- 
strated for the Be atom and then also for a few other 
atoms and small molecules [l^ . 

In the previous paper [l7| we have outlined a proof 
showing the two nodal domains property for ID har- 
monic fermions in a closed-shell singlet state for arbitrary 
size. The most interesting feature of the proof was that 
the result was very robust in the sense that almost any 
arbitrary weak interaction would induce the change in 
the topology of the nodal surfaces. 

We refresh some of the key notions and arguments 
here. Let us assume a spin-unpolarized system in its 
closed-shell singlet ground state of IN particles. Con- 
sider a simultaneous exchange of an odd number of spin- 
up pair(s) and an odd number of spin-down pair(s). For 
noninteracting wave functions such simultaneous pair ex- 
changes imply that the node will be crossed once or mul- 
tiple times. This must be the case whenever the spin-up 
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and -down subspaces are independent of each other, such 
as in the mean-field and HF wave functions. However, if 
there exists a point Rf such that during the simultaneous 
spin-up and -down pair exchanges the inequahty |^| > 
holds along the whole path, then the wave function has 
only two nodal cells. 

Using this property we can now demonstrate the two 
nodal cell for several types of systems. Under rather 
general conditions, we will show that the correlation in- 
cluded in the Bardeen-Cooper-Schrieffer (BCS) pairing 
wave function 0, ^3 given by 

*BCs(l,...,2A^) =det[$(z,j)] (29) 

smoothes out the noninteracting four nodal cells into 
the minimal two. Here <&(i,j) = is a singlet 

pair orbital for i ] and j [ fermions and we decom- 
pose it into noninteracting and correlated components 
<i>(i,j) — $o(*, j) + 'I'corr(«, j). Using one-particlc or- 
bitals we can write $o(*jj) = '^a'^a{'>')4'a{i) where 
the sum is over HF (or noninteracting) orbitals while 

^corrihj) = J2af3'^a(34>a{i)(t>f3{j) whcrC {Ca/3 = Cf3a} are 

variational parameters and the sum is over unoccupied 
or virtual ("correlating") orbitals. The BCS wave func- 
tion was originally introduced for conventional supercon- 
ductors, however, it proved to be very successful also in 
describing the correlation effects in electronic structure 
problems 0,E|. 

VI. a. 3-D harmonic spin-unpolarized fermions with 
interactions. 

For the 31? harmonic oscillator one can show that the 
interactions lead to the minimal number of nodal cells in 
a rather simple and elegant way. First let us consider a 
small system which is easy to analyze. We illustrate the 
idea on 2N — 8 particles in the singlet ground state with 
the particle positions given in Fig. [7| The correspond- 
ing closed-shell singlet ground state for the 3D harmonic 
oscillator is ^ S{ls'^2p^). (Note that for the harmonic po- 
tential the 2p state is below the 2s state, unlike for the 
Coulomb potential.) For simplicity, we drop the gaus- 
sian prefactors since they do not affect the nodes and the 
pairing functions can be then written as follows 

$o(j,i) - l + 3r, -r, (30) 

and 

'i>corr(«,j)=7[3(r,T,)2-r2r2] (31) 

where the noninteracting part is constructed from the oc- 
cupied Is, 2p orbitals while for the correlating part the 
unoccupied 3d subshell orbitals were used. The constant 
7 is a variational parameter. Clearly, both the pair or- 
bitals and the wave function are spherically symmetric. 



Assuming the positions as given in the Fig. [7| which are 
of high symmetry and therefore simplify the evaluation 
of the wave function, we get 

^'(1,...,8) = 47cos(p[l - (3 + 7)cos2(p]/3 (32) 

where ri = rs = 0, Vi = l,i = 2, 3, 4, 6, 7, 8 and cos ip — 
r4 • rs. 
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FIG. 7: Positions of 2N = 8 fermions of the 3-D harmonic well 
system. Fermions 1 to 4 have spins up, while 5 to 8 have spins 
down. Fermions 1 and 5 are at the origin while the rest of 
the particles lie on the surface of the sphere with r = 1. Note 
that rotation of the system by tt around the axis exchanges 
particles 3, 4 and 7, 8. See the text for details. 

For 7 = the wave function vanishes since the parti- 
cles in both spin subspaces lie on the noninteracting node 
(ie, four particles on the same plane). The key point is 
that for 7^0 and ip ^ tt/2 the wave function does not 
vanish. Since the wave function is spherically symmetric 
we can rotate the particles around the z— axis by tt. This 
transformation exchanges particles 3 and 4 in the spin-up 
channel and particles 7 and 8 in the spin-down channel. 
The wave function is rotationally invariant and nonvan- 
ishing what clearly implies that the BCS wave function 
has spin-up and -down subspaces interconnected since si- 
multaneous exchanges in the spin-up and spin-down does 
not hit the node. On the other hand, it is easy to check 
that the rotation of the particles in one spin channel only 
causes a node-crossing since then cos vanishes either at 
7r/2 or 37r/2. 

It is clear that our argument is correct whenever the 
strength of the interaction is small so that the BCS wave 
function is accurate enough. The variational parameter 
7 is related to the interaction strength. This implies that 
the topological change from two to four nodal cells takes 
place for arbitrary small, but nonvanishing, interaction 
strength. The only assumptions were that the interac- 
tion will induce electron correlation and lead to a nonzero 
variational parameter 7 (or, in general, to a nonvanish- 
ing correlation component) and that the wave function is 
spherically symmetric. Our analytic argument therefore 
shows that for weak interactions the nodal cell "degener- 
acy" is lifted and the multiple cells are smoothed to the 
minimal number of two. 
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It is important to note that one can find difi^erent pic- 
tures in other circumstances. For example, for strong 
or nonlocal interactions, imposed additional symmetries 
or large degeneracies, the nodal changes can be differ- 
ent and the resulting nodal topology might exhibit more 
than the two nodal cells; this aspect is further discussed 
in the conclusion section and require more investigation. 

Coming back to our example of harmonic fermions with 
interactions, using similar arguments the two nodal cells 
can be demonstrated for an arbitrary size closed-shell 
ground state. Consider 2N particles in a singlet closed- 
shell ground state. The closed-shell states can be labeled 
by M which represents the "Fermi momentum" for the 
harmonic oscillator and N = {M +1){M + 2){M + 3)/6. 
Let us decompose N into the maximum odd number of 
pairs Np and the rest. Therefore we write N = 2Np + J 
where Np is an odd integer while J has one of the values 
J = 0,1,2,3, so that we can also define J — NmodA. 
For example, in the previous example with 2N — 8 par- 
ticles we have Np — 1 and J — 2. First, we place J 
particles from each spin channel on the z-axis in distinct 
positions. Next, we form Np pairs in each spin sub- 
space so that, say, the pair i |, (i -I- 1) | is positioned as 
given by {ri,'di,ipi) = {rk,i}k, ^Pk), {n+i^-di+i, ipi+i) 
(''fc, ^fc, <y5fe -t- tt), where r^d,ip are the spherical coordi- 
nates. Here fc = 1, ...jNp labels the pairs in the spin-up 
channel; the spin-down particles are placed similarly and 
labeled by the pair index / — Np + 1, ...,2Np. In this 
configuration the particles lie on the noninteracting node 
since 



det [$o(«,j)] = det 



<i)nm{i)(t>nm{j) 



n+m<M 



det [(l)nm{i)] det [(l)nm{j)] = *HF*HF 



(33) 



The rotation of the system by tt causes node-crossings 
in both spin channels so that both spin-up and -down 
Slater determinants must vanish due to the rotational 
invariance. Now, if all the Np pair distances and angles 
rk,ipk7i!^k,'r'i,'di,Lpi are distinct, then each of the matri- 
ces [0nm(*)] , [<t>nmU)] ^as cxactly One linearly dependent 
row, ie, both have ranks Np — 1. This can be veri- 
fied directly for small values of M and then using in- 
duction for any M. Consequently, the matrix [$o(*jj)] 
has linear dependence in one row and one column, ie, 
it has the rank of Np — 1 as well. In general, adding 
virtual states through ^corrihj) provides independent 
rows/columns which eliminate the linear dependency so 
that det [^o{i,j) + ^corrih j)] is nonzero. Let us now as- 
sume that the interactions do not break the rotation in- 
variance. Since at this point the wave function does not 
vanish and is rotationally invariant, it does not vanish 
for the whole exchange path implying that the correlated 
BCS wave function has only two nodal cells regardless of 
the size 2N. 



Clearly, the assumption of the rotational invariance for 
the interacting case might be too restrictive since the in- 
teractions could possibly break the spherical symmetry. 
To demonstrate the two nodal cells in such a case one 
would need to show that the wave function does not van- 
ish for the complete exchange path. 



VI. b. 2D harmonic spin-unpolarized fermions with 
interactions for A'^ — 2Np + J when J = 2, 3. 

In our previous paper we have demonstrated the two 
nodal cells for the ground states with the number of par- 
ticles 2N where N = 2Np + J and J = 0, 1. (We use the 
same notation as in 3D where Np is the maximum odd 
number of pairs and J = 0, 1, 2 or 3). Note that the states 
with J = or 1 exist at any size. Here we would like to 
present an alternative and a little bit longer proof which 
covers J = 0, 1 but also remaining cases when J — 2,3, 
eg, for iV = 21 = 2 X 9 3, iV = 28 = 2 x 13 -I- 2, 
etc. In 31), the rotation axis can accommodate multiple 
particles so that one can always form an odd number of 
pairs for the exchanges. In 2D, the symmetry point (ori- 
gin) can accommodate only one particle from each spin 
subspace therefore the proof has to be modified. For 
this purpose we outline a factorization which is different 
from the previous paper |l7| . Let us remind that for 2D 
harmonic oscillator the closed-shell states and the sys- 
tem size are labeled by M = 1, 2, ... where n + m < M 
with the number of fermions in one spin channel given 
by = (M -I- 1)(M -I- 2)/2. We express the one-particle 
states as polynomials in and (re**^)™ where r, (p are the 
cylindric coordinates (we again omit the gaussian factors 
which are irrelevant for the nodes) . Note that for a given 
M the quantum number m increases or decreases by 2 
unlike in 3D. Therefore for M = 2 we have states re^"'', 
for M — 3 the additional states arc (r^ — 1), r^e'^^''', etc. 
We first write down the wave functions for three- and 
six-particle spin-polarized systems with positions given 
in Fig. IHl 

The three particle wave function is given by 

det{l,re±*'^} 

^2d{1,2,3) = ^lorl Y[ sin((^,fc/2) (34) 

l(2)<j<fc<3 

where fiQ is a constant. The product lower bound is 2 if 
the positions of the particles are as in Fig.|Sli, while it is 
1 if the positions are as in Fig. |SJd. Note that in the case 
of configuration in Fig. [SJi rotation by tt flips the wave 
function sign so that the wave function vanishes while for 
the configuration in Fig.|Ht> it does not. For six particles 
in the positions outlined in Fig. (Sji the wave function is 
clearly nonvanishing since we can write 



2D 



(1, 6) - - n sin(¥',fc/2) (35) 



3<k 
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FIG. 8; Alternative positions of three and six spin-polarized 
fermions for the 2D harmonic potentiaL (a) Three particles 
positioned on the noninteracting node, (b) Three particles in 
positions with a nonvanishing wave function, (c) Six particles 
positioned on the noninteracting node, (d) Six particles in 
positions with a nonvanishing wave function. See the text for 
details. 



where /ig is some constant (in Fig. [SJi we have ri = 0) . 
Alternatively, the six particles can be positioned in pairs 
as sketched in Fig. |SJ: where the particles are positioned 
on the HF node so that the noninteracting wave function 
vanishes. 

For the sake of completeness we provide the expression 
for the wave function for a general system with size M 



*2D(l,...,iV) 



M 

)n 

fc=i 



n 

l<j<k 



(36) 



where M = int[(A/+l) /2] is the integer part of (A/+1) /2. 
The particles with indices in the subset Ik lie on a circle 
with the radius and Uk = 2{M + 1 — 2k) + 1. The 
prefactor is again a nonnegative function of Tj which has 
no impact on the nodes and is constant for rotations. 
From the last equation we see that in cases when M is 
even, eg, M — Q,N — 21, we end up with factorization 
of the type 11 x 7 x 3 particles while for M odd, eg, 
M = 7,7V = 28, we get 13x9x5x1. That means that the 
last three particles (M even) or the last six particles (Af 
odd) can be always arranged into possibilities as sketched 
in Fig|SJi-d. It is now clear that we can prove that the 



spin-up and -down configuration space are interconnected 
for the cases when J = 2,3 as specified above. Using 
the alternative configurations for three and six particles 
we can arrange the systems into such positions that the 
rotation by tt exchanges odd number of pairs in both 
spin channels. For example, if N is odd and [N — l)/2 
is even we position three particles on a triangle as given 
in Fig.lHt). That takes out one pair from the (A^ — l)/2 
exchanges and the proof then follows using the arguments 
for ZD harmonic fermions. Similarly, if TV is even and 
A^/2 is also even, we position one particle of each spin 
at the origin and five from each spin channel on a circle. 
This eliminates three pairs from 7V/2 and again we end 
up with remaining number of pairs being odd so that the 
rest of the proof follows similarly to the ZD case. 



VI. c. Spin unpolarized interacting 2D and 3Z) 
homogeneous electron gas. 

Let us now prove the two nodal cells for the spin- 
unpolarized closed shells for d > 1 homogeneous gas 
with interactions. We use the relevant definitions intro- 
duced previously for the spin-polarized periodic fermion 
gas (section III. a.). We assume a 2D system of 2N par- 
ticles in the spin singlet ground state where the one- 
particle states occupy the Bloch states within the Fermi 
disk with the radius kp- Let us first consider 2iV = 10 
particle case and we position the particles as in Fig. 
The wave function is translationally invariant therefore 
the action of leaves the wave function unchanged. As- 
suming first that there is no interaction, we see that the 
translation exchanges the particles 1, 2 but due to the in- 
variance the wave function value does not change: that 
is possible only if the particles are sitting on the node. 
This also agrees with expressions for the wave functions 
derived before (see Eq. 0) Consider now that we switch- 
on interactions and describe the particles by the BCS 
wave function as given above. The pairing orbital for the 
noninteracting gas is given by 

*o(«,i)== X! 0k(i)0-k(j) = l + 2[cos(xij) + cos(?;y)] 

|k|<feF 

(37) 

where we used the occupied orbitals {1, e^'^, e^'^}. The 
correlated part is given by 

$corr(i, j) = 7 X! '^k(i)0-k(j) = 
|k|>feF 



7[cos(xy ) cos(yy ) -t- coii{2xi.j) + cos{2yij)] 



(38) 



where ^corr{hj) includes the next two unoccupied stars 
of states with k = v^, 2 and 7 is a variational parameter. 

Unlike for the noninteracting wave function (ie, 7 = 0) 
the BCS correlated wave function does not vanish and 
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FIG. 9: Positions of five spin-up electrons from an interact- 
ing ten-fermion system in the 2D periodic box (— 7r,7r)^. For 
clarity, the five down-spin particles are not shown. The spin- 
down particles are positioned in the same pattern on two lines 
X — const. 



this is indeed the case due to the arguments presented 
above. The Slater matrix of uncorrelated wave function 
has linear dependency while for the correlated case there 
always exists a configuration of lines with nonvanishing 
wave function. This is due to the elimination of the lin- 
ear dependency through addition of virtual orbitals, as 
explained for the 31? harmonic oscillator. Note that as 
soon as the wave function is translationally invariant the 
wave function does not vanish for the whole translation 
path, implying that the spin-up and -down nodal do- 
mains are interconnected and the wave function has only 
two nodal cells. It is straightforward to extend the proof 
to arbitrary size closed-shell system. The configuration of 
particles can be given as follows: an odd number of pairs 
(eg, one) in each spin channel is positioned in a similar 
way as in so that the translation by tt exchanges the 
particles in the pair. For example, for a spin-up pair 
i,i + 1 we specify the coordinates as yi = — j/i+i = 7r/2 
and Xi — Xi-i-i — ^0 where ^0 is distinct from x coordi- 
nates of the rest of particles with the same spin. The rest 
of the particles is positioned in such a way so that the 
brings them into a symmetric position regarding the 
reflection around axis y. The wave function is transla- 
tionally invariant and therefore uncorrelated wave func- 
tion vanishes while for BCS case it is, in general, nonzero. 
The two nodal cells property then follows using the same 
general arguments as in the previous cases. 

The correlated wave function has another important 
property namely that one can wind around the periodic 
box a singlet pair of particles (ie, spin-up and spin-down 
pair) without hitting the node. Indeed, this is implied 
by the fact that there are only two nodal cells: once this 
is the case one can then always find such a path that 
the pair of particles can wind around the box without 
node-crossing. For example, for the ten-particle exam- 
ple above we can wind spin-up and -down particle pair 
around along y = Co without hitting the node. We plot 
the wave function for two types of particle positions: first, 
the spin-down particles are at identical positions as the 



spin-down ones, ie, rg = ri, ry = r2, etc fFig. llOb): sec- 
ond, the particles in spin-up and -down channel were off- 
set by small displacement ri+5 = r.^ + 0.2, i 1, 5, see 
Fig.llOb. Consider now the simultaneous translation r^*^ 
of the particles 1 and 6. This translation winds the pair 
around the box and the wave function values for the paths 
are plotted in Fig. llOb.b. The value of the displacement 
between the particle positions neither the repositioning 
of the particles around the box is not crucial and there 
exists a significant size subspace of particle positions for 
which the winding of the pair is possible. We have chosen 
just two simple examples as qualitative illustrations. 

Note that the wave function enables to wind also a sin- 
gle particle without the node crossing. However, single 
particle winding is possible both for correlated and un- 
correlated wave functions since, say, the spin-up HF part 
has two nodal domains and correlated wave function only 
smoothes out the HF nodes. Therefore this property is 
not affected by the correlation in any significant manner. 

The effect that singlet pair of particles can pass 
through nodal openings between the spin-up and -down 
subspaces have been demonstrated for small number of 
particles before, for example, in our paper on employing 
pfaffi an p airing functions for electronic structure prob- 
lems [lal- The case here illustrates that this property 
stems from interconnected spin-up and -down subspaces 
and therefore applies to similar two nodal cell wave func- 
tions, in general. 

It is well known that a Fermi liquid, such as the ho- 
mogeneous electron gas, becomes unstable to a weak at- 
tractive interaction, develops Cooper pair instability and 
opens the superconductivity gap [l2|. This effect is cap- 
tured by the BCS wave function which we used in the 
proof above. Cooper pairs can therefore wind around the 
box without hitting the node although this effect is not 
exclusive to Cooper pairs only. In fact, the connectedness 
of spin-up and -down subspaces is a rather generic prop- 
erty in the sense that it appears also in systems which are 
not necessarily superconducting, ie, in Fermi liquids with 
repulsive interactions. Since superconductivity is charac- 
terized by macroscopic phase coherence and a number of 
other properties, which might or might not be related 
to the nodal topologies, on the basis of the analysis of 
the BCS wave function above one expects that the nodal 
opening necessarily appears in the superconducting state, 
however, for the superconductivity to occur this condi- 
tion is not sufficient. In addition, here we assume only 
the simplest s— wave pairing; for the wave or higher 
angular momentum pairing the situation is less clear and 
needs to be further investigated. 
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FIG. 10: Winding the particles 1 (spin-up) and 6 (spin-down) 
around the 2D box along the hue y — (^o- The dashed hne 
is the Hartree-Fock wave function (arb.u.) while the full line 
is the correlated BCS wave function, a) Particles 1 and 6, 2 
and 7, etc, share the same position. The uncorrelated wave 
function touches the node quadratically since the spin-up and 
-down determinants are identical, b) The positions of parti- 
cles in spin-down channel are offset by 0.2 from the positions 
in the particle spin-up channel. The uncorrelated wave func- 
tion crosses the node multiple times since the spin-up and 
-down determinants are different. 



VII. NODES OF FERMIONIC DENSITY 
MATRICES. 

In this part we generalize the ideas presented in pre- 
ceding section to temperature dependent density matri- 
ces |5j. Consider first a system of spin-polarized fcrmions. 
The temperature/imaginary time density matrix is given 

by 

g{R', R,f3)=J2 e-'^^"^:(i?')*n(^) (39) 

n 

where /? is the inverse temperature and the sum is over 
the complete system of eigenstates of a given Hamilto- 
nian H. It is clear that the density matrix is antisym- 
metric in particle exchanges in the same manner as the 



wave function ^{R) or ^{R'). Therefore the notion of 
fermion nodes can be generalized also to density matrix 
in {2dN + 1) dimensions since there is an explicit de- 
pendence on P as well. As pointed out elsewhere |3, for 
fixed R' and /3 one can study the nodes and nodal cells in 
the dN dimensional i?— subspace. Similarly to the wave 
function, the node then becomes {dN — l)-dimcnsional 
manifold with the generalization that it is dependent on 
R' and (3. For the fixed R' and (3 the tiling property 
holds in the same manner as for the wave functions. The 
key additional feature of the density matrix is that once 
there are only two nodal cells at some initial f3o than this 
property holds for any (3 > . This is not difficult to 
understand since the density matrix fulfills the following 
linear equation 

dg(R,R',P) , , , 

- Qp =Hq{R,R',P) (40) 

with an initial condition 

q{R, R', 0) = Ad{R - R') = det[(5(r, - r'^)] (41) 

where A is the antisymmetrizing operator. 

We now understand that for almost any Hamiltonian 
with interactions the density matrix will have only two 
nodal cells for sufficiently large /3 (ie, at low tempera- 
tures). This is due to the fact that at sufficiently low tem- 
perature the ground state becomes dominant (Eq. I39|l . 
The key point now is to show that this is the case also 
for high temperatures. The free particle density matrix 
is given by 

giR,R',P) = (27r/3)-'^^/2det [exp(~|r, - r',f /2/3)] 

(42) 

where we assume atomic units with h^/m — 1. (For 
other than free boundary conditions, such as for the pe- 
riodic ones, the expression has to be modified accord- 
ingly.) Note that this density matrix is universal since at 
sufficiently high temperatures the interactions become ir- 
relevant. 

There are several ways how to prove that the high- 
temperature density matrix has only two nodal cells. For 
very small /? one can use the induction as follows. As- 
sume that TV particles are described by the density matrix 
given by Ea.E^and the particles are connected by triple 
exchanges for a fixed i?',/3. We add an additional par- 
ticle with the label iV -I- 1 to the system which occupies 
a certain region of the configuration (free) space. The 
particle TV -I- 1 is positioned at the border of the occupied 
region and let us assume that the particles with labels 

— 1 and N are its closest neighbours. Let us move 
the three particles N — 1,N,N + 1 away from the rest 
without crossing the node (what can be always done by 
appropriate positioning). Since for small /3 the overlaps 
of gaussians become small, one can factorize the determi- 
nant into a product of the three particles N — 1, N, N +1 
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determinant and the determinant for the rest. For the 
three particles the density matrix has only two nodal 
cells as one can show easily and therefore the addi- 
tional particle is connected. This applies to both R and 
R' subspaces since they must have identical properties. 
(In what follows we will show that, in fact, at high tem- 
peratures the primed and unprimed spaces are connected 
as well.) 

There is also an alternative proof which is interesting 
also on its own since it provides a different view on the 
density matrices. Note that the functional form of the 
free particle density matrix has a unique property. Com- 
paring the Eq. 021 with the BCS wave function (Eq. 1221) 
we see that the high-temperature density matrix can be 
identified with a BCS wave function if we properly de- 
fine an underlying effective model. Instead of our origi- 
nal system of N spin-polarized fermions, let as consider 
a model system with 2N particles so that the configura- 
tions R and R' denote positions of these different sets of 
particles, which we will call for simplicity unprimed and 
primed particles. Let us then define a new Hamiltonian 
H(R, R') with an effective quadratic interaction between 
the unprimed and primed particles as given by 



HiR, R') = T + T' + Vof3-' ^ |r. - r'. 



(43) 



"excited states" (ie, beyond the Fermi sphere) due to the 
primed-unprimed interactions, we find that the unprimed 
and primed nodal cells are interconnected. Therefore at 
classical temperature Po the density matrix has only two 
nodal cells and then the same is true for arbitrary P > (3q. 
This primed-unprimed interconnection becomes less pro- 
nounced and ceases to exist at /? — > cxd since then the 
density matrix is proportional to the " noninteracting" 
product ^'o(i?)^o(-R') where "^o is the ground state for 
the original physical system of N fermions described by 
H. This proof is therefore based on an interesting duality 
between the A'^ spin-polarized fermions at classical tem- 
peratures and the model system with 2N particles with 
a temperature-dependent, harmonic interaction between 
the unprimed and primed subspace particles. 

Possibly, there might be also a third way how to 
prove the minimal number of nodal cells for the high- 
temperature density matrix through diagonalization of 
the quadratic Hamiltonian H (we have not investigated 
this possibility). The fact that also the density matrices 
have two nodal cells is important for the path integral 
Monte Carlo methods which for fermions often employ 
the fixed-node approximation adapted for the path inte- 
grals 113. 



where T and T' denote kinetic energy operators for the 
corresponding sets of particles and Vq is a constant with 
appropriate dimensions. Note that particles within the 
given set, say, unprimed, are antisymmetric but other- 
wise do not interact with each other. The interaction 
appears only between the primed and unprimed degrees 
of freedom as given by the Hamiltonian H{R,R'). For 
(3^0 the exact wave function for this system is a BCS 
wave function \E'(i?, R') = det[(p{ri, r'j)] where the indices 
i and j label unprimed and primed particles, respectively. 
The pairing function (f)(Ti,r'j) is obviously the gaussian 
given above. It is not too difficult to demonstrate that 
this wave function (and the density matrix) has only two 
nodal cells. For example, we can expand the gaussian 
into plane waves 



exp(- 



V2/3) 



k 



(44) 



where {ck} are expansion coefficients. The sum is over 
states within the Fermi sphere of a periodic box which 
accommodates 2N particles for a given density. Let us 
specify that N is large and (3 is such that the system 
can be considered classical so that the actual interaction 
in the original Hamiltonian H is irrelevant. Then the 
density matrix corresponds to the Hartree-Fock product 



det[e*'^'-'"^]det[e-''''"''^] 



(45) 



However, we have already proved that such system has 
only two nodal cells. In addition, if the sum includes also 



VIII. DISCUSSION AND CONCLUSIONS. 

Inspired byprevious conjectures and numerical studies 
the presented analysis and proofs gen- 
eralize and clarify the properties of ground states fermion 
nodes. We have employed symmetries, wave function fac- 
torizations and triple exchanges to prove that for d > 1 
the closed-shell ground states have two nodal cells for ar- 
bitrary number of particles in several paradigmatic mod- 
els. In this paper the proofs were carried for closed shells, 
mainly to avoid additional complications from degenera- 
cies and some of these aspects will be addressed in sub- 
sequent papers. 

It is perhaps more interesting to discuss viable possi- 
bilities when the ground state wave functions might have 
more than two nodal cells. Clearly, it is easy to gener- 
ate more nodal cells by imposing additional symmetry 
or boundary conditions. Another possibility comes from 
very strong interactions, for example, whenever the ef- 
fect of interaction would become competitive with the 
kinetic energy increase from forming additional nodes it 
is possible that more than the minimal two nodal cells can 
form. One can also construct nonlocal interactions which 
would violate some of the properties mentioned here (eg, 
the tiling property), or reorder the states so that, for ex- 
ample, excited states could lie below the ground state, 
etc. Another candidate for unusual effects in the nodal 
structure are open-shell systems with large degeneracies 
and densities of states at the Fermi level. Some of these 
interesting issues will be subject of future studies. 
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In conclusion, we have presented a number of new re- 
sults which reveal the structure of nodes and nodal cells 
of fermionic wave functions. Building upon ideas intro- 
duced in previous paper we were able to demonstrate 
the minimal number of two nodal cells in several spin- 
polarized models such as noninteracting fermions in a 
periodic box, in a box with zero boundary conditions, 
fermions on a spherical surface, etc. This enabled us to 
formulate a theorem which states that in d > 1 the min- 
imal two nodal cells are present for any Slater determi- 
nant with monomial matrix elements for any size which 
allows for a closed-shell nondegenerate ground state. We 
have shown that this property extends also to cases which 
cannot be described by the Slater matrix of monomials 
such as the noninteracting and HF atomic states up to 
the 3d shell. We have studied the effect of interactions on 
the noninteracting nodal cells of spin-unpolarized, closed- 
shell singlets. Our results show that, in general, the in- 
teractions smooth out the noninteracting multiple nodal 
cells into the minimal number of two. For interacting ho- 
mogeneous electron gas we have demonstrated that the 
two nodal cells allow singlet pairs of particles to wind 
around the periodic box without crossing the node. Fi- 
nally, we have shown that the temperature/imaginary 
time density matrix has very similar nodal structure and 
therefore the minimal two nodal cells property applies 
also to density matrices with important implications for 
path integral Monte Carlo simulations. We have demon- 
strated this by using an appropriate mapping of the den- 
sity matrix onto the ground state of a model systems with 
twice as many particles interacting with temeperature- 
dependent harmonic potentials. 
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